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Abstract 

By employing the methods of renormalized field theory we show that the 
percolation behavior of random resistor-diode networks near the multicritical 
line belongs to the universality class of isotropic percolation. We construct a 
mesoscopic model from the general epidemic process by including a relevant 
isotropy-breaking perturbation. We present a two-loop calculation of the 
crossover exponent <f>. Upon blending the e-expansion result with the exact 
value 4> = 1 for one dimension by a rational approximation, we obtain for two 
dimensions <j> = 1.29 ±0.05. This value is in agreement with the recent simula- 
tions of a two-dimensional random diode network by Inui, Kakuno, Tretyakov, 
Komatsu, and Kameoka, who found an order parameter exponent (3 different 
from those of isotropic and directed percolation. Furthermore, we reconsider 
the theory of the full crossover from isotropic to directed percolation by Frey, 

Tauber, and Schwabl and clear up some minor shortcomings. 
PACS-numbers: 64.60.Ak, 05.40.-a, 64.60.Ht, 64.60.Kw 
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I. INTRODUCTION 



Random resistor-diode networks (RDN) were introduced by Redner []!]-§]. Neverthe- 
less, they were already contained implicitly in the pioneering work of Broadbent and Ham- 
mersley [|J on percolation. RDN define a percolation model (for a recent introduction to 
percolation see Stauffer and Aharony [||) on a d- dimensional hypercubic lattice in which 
nearest-neighbor sites are connected by a resistor, a positive diode (conducting only in a 
preferred direction), a negativ diode (conducting only opposite to the preferred direction), 
or an insulator with respective probabilities p, p + , and q = 1 — p — p + — p_. In the 
three dimensional phase diagram (pictured as a tetrahedron spanned by the four proba- 
bilities) one finds a nonpercolating and three percolating phases. The percolating phases 
are isotropic, positively directed, or negatively directed. Between the phases there are 
surfaces of continuous transitions. All four phases meet along a multicritical line, where 
< r := p + = p_ < 1/2 and p = p c (r). On the entire multicritical line, i.e., independently 
of r, one finds the scaling properties of usual isotropic percolation (r = 0). 

About 20 years ago Redner [|l]-|3|] studied the phase diagram sketched in Fig. 1 as well 
as geometrical properties of RDN in two dimensions. He used real-space renormalization 
methods and planar lattice duality ||. 

Recently Inui et al. [|7| have measured the order parameter exponent (3 for the special 
case of a two-dimensional random diode network (with p = q = 0) by Monte Carlo methods 
combined with series expansions. At the symmetric critical point p + = p_ = 1/2 (so called 
random Manhattan (RM)), they found j3 = 0.1794 ± 0.008, which does not coincide with 
the known values wether for isotropic percolation (IP), (3ip = 5/36, nor directed percolation 
(DP), f3r)p = 0.27643. Therefore, they concluded that the percolation properties of the ran- 
dom diode network constitute a new universality class different from isotropic and directed 
percolation. 

In this paper we study RDN by the methods of renormalized field theory. Contrary to 
Inui et al. we find that random diode networks at the percolation point as well as RDN at the 
full multicritical line belong to the universality class of isotropic percolation. The variable 
r = p + = p„, which maps out the multicritical line p = p c (r), is a redundant variable in the 
sense of the classification scheme of scaling variables by Wegner ||. Thus, all the points of 
this line are equivalent to the usual isotropic percolation point with r = 0. 

The case r > 0, in which positive and negative diodes are distributed with equal proba- 
bility, leads to a breaking of isotropy and to elongated percolating clusters. However, this 
symmetry breaking can be easily compensated in the mesoscopic field theoretic formula- 
tion by a simple rescaling of the length scale of the preferred direction. The specific length 
scale is therefore redundant. Wegner showed that the renormalization flow of a redundant 
variable depends on the particular form of renormalization group used and does not affect 
the physics. Thus, in the case of Redner's real space renormalization group, a special fixed 
point, the so called mixed one, is distinguished. 

The situation is different for a symmetry breaking which favors not only an axis but also 
a direction on that axis. This leads to a relevant variable ~ (p + — p_) and therefore to a new 
correlation length exponent v = (pi^ip. The crossover exponent describes the beginning 
of the crossover to directed percolation. It also relates the order parameter exponents, 
P = (f>Pip. In two dimensions, it has, according to the results of Inui et al. J7|, the value 
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= 1.29... . Our perturbation calculation yields <p = 1-29 ± 0.05 in agreement with the 
simulations. It should be noted that the interpretation of the multicritical line as a line 
of equivalent fixed points with the same scaling behavior, as well as the new correlation 
exponent u, can be already found implicitly and explicitly in the early papers of Redner 
]l|-[3[ (his notation is v + for the new exponent v). 

The organization of the remainder of this paper is as follows. In Sec. II, we first develop 
a mesoscopic field theoretic model that is capable to describe the crossover from isotropic to 
directed percolation, which is the basic feature of RDN. We describe the renormalization of 
the model and calculate the renormalization factors to two-loop order. In Sec. Ill, we intro- 
duce the renormalization group equation for the model and derive the general asymptotic 
scaling properties. In Sec. IV, we derive an interpolating formula for the crossover exponent 
from the two-loop ^-expansion and show that it reproduces the result of Inui et at In Sec. V 
we reconsider the theory by Frey, Tauber, and Schwabl [[| for the crossover from isotropic 
to directed percolation and clear up some minor shortcomings. In Sec. VI, we give some 
conclusions and summarize our work. In an appendix, we present the two-loop calculation 
of the renormalization factors. 



II. THE FIELD THEORETIC MODEL AND ITS RENORMALIZATION 



Here we develop a mesoscopic model that is capable to describe the crossover from 
isotropic to directed percolation. In this paper we are only interested in connectivity prop- 
erties of the percolating system. In contrast to earlier work on random resistor networks 
T0|-|T3| we neglect transport properties as the conductance etc. In other words: all we ask 



is wether two points on the lattice are connected or not. Formally, we consider the limit of 
zero resistance of the conducting elements. 

Percolating clusters in space and time can be generated by a stochastic spreading process 
known as the general epidemic process (GEP) |JL4|| . In order to apply field-theoretic methods 
it is convenient to use the path-integral representation of the underlying stochastic 
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process s(x, t) [|I7HT9"11. With the imaginary- valued response field denoted by s(x, t), the gen- 
erating functional of the Greens functions, the connected response, and correlation functions 
takes the form 



W[H, H} = ln V{s, s) exp 



J[s, s}+ d d x / dt(Hs + Hs 



(2.1) 

The dynamic functional J~\s, s] and the functional measure T>(s, s) (T>(s, s) is a symbolic 
notation for Y\ xt (ds('x,t)ds('x,t)) times a constant) are understood to be defined using a 
prepoint (Ito) discretization with respect to time [|19|1 . The prepoint discretization leads 
to the causality rule 9(t < 0) = in response functions. This causality rule then forbids 
response propagator loops in the diagrammatical perturbation expansion. 

Using the path-integral formulation a renormalized field-theory of dynamic isotropic 
percolation can be gained from the GEP [^(],[np. If we add to this model a relevant coupling 
which breaks isotropy and introduce a preferred direction n for the spreading of the disease 
in the <i-dimensional space, we get the dynamic functional 



J 



d xl / dts 



d t + \(t - V 



2v(n ■ V)) - ^-s 



s H so t 



dt's(t') 



(2.2) 



3 



This functional corresponds to the Langevin equations 

^s(x, t) = A(V 2 - 2v(n • V) - r - <?n(x, i))s(x, t) + C(x, t) , (2.3) 
n(x,t) = A f s(x,t')dt' , (2.4) 

J — oo 

(C(x, t)C(x', 0) = A^(x, t)5(x - x')<5(t - , (2.5) 

for the GEP with the suitable scaled density s(x, £) of the infected individuals. Here n(x, t) 
constitutes the density of the immun (or dead) individuals and C( x )^0 is a Gaussian noise 
which is zero in spatial regions where the disease is extinguished. The deterministic drift of 
the disease in space is represented by the flow 2\v ns(x, t). 

Of course, if the d- dimensional rotational symmetry is broken to (d — 1) -dimensional 
isotropy, the diffusion constants for longitudinal (with respect to the direction n) and 
transversal spreading are in general different. Thus we have to consider a more general 
diffusion operator AV 2 — > A(V^ + c~ 2 <9jj). However, it is easy to see that the new param- 
eter c can be absorbed into the definition of the longitudinal length scale by cx\\ — > x\\ 
followed by an appropriate change of the densities and the coupling constant g. After that, 
diffusional spreading looks isotropic again. The parameter c depends on the microscopic 
model. It is a redundant variable in the sense of Wegner H and is responsible for the mul- 
ticritical line in the RDN fTj-Q. From the microscopic RDN- standpoint, the three variables 
r, c, v are analytical functions of the three probabilities p, p + , p_ for the resistors and 
diodes near the critical manifolds of RDN and share their spatial symmetries. Thus one has 
t(p,P+,P-) = r(p,p-,p+), c(p,p + ,p_) = c(p,p-,p+), and u(p,p + ,p_) = -v(p,p-,p+). In 
particular we have c(p, 0, 0) = 1, but in general c(p, r, r) ^ 1 if r > 0. Moreover, v(p, r, r) = 
holds. Remember that p + = p_ = r and p = p c (r) defines the multicritical line. 

In this paper we are interested only in the static behavior of the process, i.e., in the 
statistics of the distributions of immunes in space n(x, oo) after a long time when the 
epidemic is extinguished. These distributions constitute isotropic percolating clusters at 
the critical point of the GEP, which is given by v — v(p c (r), r, r) = and r = r(p c (r), r, r) = 
r c = if we neglect fluctuation corrections. The density n(x, oo) is proportional to the 
Fourier transform of s for frequency zero. The statistical weight for these frequency-zero 
modes can be found by invoking the formal limit 

/oo 
dts(x,t)=n(x, oo) —>(p(x) , (2.6) 
-oo 



for the fields in the dynamic functional J~, Eq. ( |2.2| ). This manipulation can be controlled 
term by term in a diagrammatical perturbation expansion and leads to the wanted represen- 
tation of the zero frequency Green's functions as functional integrals with a weight exp(— Ti.) 
and a quasistatic Hamiltonian 



H = J d d x\(p r- V 2 + 2w(n- V) + |( 



<p-<p 



V ■ (2-7) 



Remember that closed loops of the response propagators are forbidden in the Feynman dia- 
grams as a consequence of the causality rule. Therefore, the functional TC is, up to a rescaling, 
identical with the statistical functional considered by Frey, Tauber, and Schwabl p|,|22 
their work on the crossover from isotropic to directed percolation. 
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To absorb ultraviolet divergencies in a perturbational calculation of the Green's functions 
with the Hamiltonian 7i we use, in the case v = 0, the following renormalization scheme [20] : 



(p -> (p = Z 1/2 fi , ip^ip = Z 1/2 <p , (2.8) 
r^f = Z~ x Z t t\? + r c , g 2 ^g 2 = G: l Z-*Z u uif . (2.9) 

Here e = 6 - d, G £ = T(l + e/2)/(47r) d / 2 , and \x is the usual external momentum scale, 
which makes the renormalized coupling constant u dimensionless. Note that the fields (p 
and p are renormalized by the same Z-factor as a consequence of the reflection symmetry 
ip(x±,x\\) «-> — <p(x_i_, of H, which leads eventually to the equality (3' = (3 between the 
exponents characterizing the particle density and the percolation probability. The renor- 
malizations are known from percolation field theory up to three-loop order |23j. Using di- 



mensional regularization and minimal subtraction (minimal renormalization) together with 
the ^-expansion one finds f c = and to two-loop order 

^ 1 + s + ( 11 -^)(s) 2+0( « 3 )- < 2 ' 10 ' 

Zr = 1 + H + ( 9 _«)(-) , + 0(ll . ) , „n, 

Au f 59 \ fu\ 2 . . 

Z u = 1 + — + (15 - —ej + 0(u 3 ) . (2.12) 

If v ^ 0, further renormalizations are needed. Because v ~ fi, this relevant parameter 
has positive naive dimension like r. Hence, we consider it as a soft variable (which means 
that the renormalization constants do not depend on r and v) as long as v is finite and 
complete the renormalization scheme by 

r -> f = Z- 1 [z T r + Y w v 2 \ii 2 , v -> v = Z~ l Z v vyL . (2.13) 

Here we anticipate that r and v 2 get mixed under renormalization. 

To calculate Z v und Y VT , we need the (unrenormalized) propagator G(q) = (y? q y?_ q )o rimc ^. 
We define the spatial Fourier transform by <£>(x) = J_ <£> q exp(iq • x) (with the abbreviation 
J q = (27r)- d /A)andget 

G(q) = =— = * , (2.14) 

V ; r + g 2 + 2iv • q (r + t; 2 ) + (q + iv) 2 V 7 

where we used the notation v =vn. The one-loop contribution to the vertex function r^i 
(the amputated one-particle irreducible Green's function with one s- and one s-leg) is given 
by 

,2 



(q) = y I G(k + q/2)G(-k + q/2) 



Y y ds l ds 2 J^exp 



-(s 1 + s 2 )(t + v 2 ) 
-si(k + (q/2 + zv)) 2 - s 2 (k - (q/2 + w)) 5 
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(2-e)e {T + Q) 4 VV r + Q 

\{Z — e)e e 



4 

where we have defined the functions 



-?^ 2) (^)+0(g 3 )) , (2.15) 



Kf\p) = / dx^ tJ rm (2-16) 

and the abbreviation 

Q = g 2 /4 + ivq. (2.17) 
The expansion of the yet unrenormalized one-loop selfenergy E^'(q) in v 2 and £ yields 

S (1 Hq) = -^^((l + |)r + y + fv-q+iq 2 + 0(. 2 )) . (2.18) 

Using the renormalization scheme Eqs. ( |2.8| . |2.9| . [2.13| ), the renormalized vertex function I\ ;1 
is found to first order in the renormalized coupling constant u as 

u / v 2 \\ o u\ / u 



r lfl = (Z T r + Y^ 2 - ^(r + ^-) J/. 2 + 2i^„ - ^j/iv ■ q + [Z - ^-Jq 2 + 0(e°) , 

(2.19) 

from which the new renormalizations can be gathered to 0(u) as 

: ) 2 + 0(u 3 ), (2.20) 



u 




(23- 


73 \ 




+ 


i2 £ J 


' 3e 








30 


35 




H 













Y - = Te + { 3 °-'T-){i- C ) +0( " 3) - (2 - 21) 
Here we have included also the two-loop result calculated in the appendix. 



III. RENORMALIZATION GROUP EQUATION AND SCALING BEHAVIOR 

Next we explore the scaling properties of percolation in the RDN system. Scaling prop- 
erties describe how physical quantities will transform under a change of length scales. By 
the renormalization, we have introduced the arbitrary mesoscopic length scale fjT 1 . The 
freedom to choose /i, keeping the unrenormalized fields and bare parameters {f, v,g}, and, 
in cutoff regularization, the momentum cutoff A fixed, can be used to derive in a routine 
fashion the renormalization group (RG) equation for the connected Green's functions 

/ N N+N , (conn) 

G ^(w)= n^w n v^)) • (3.1) 

\=l j=N+l 1 
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We denote /^-derivatives at fixed bare parameters by 9 M | Q . From \i d^G^^^ 

renormalization scheme, Eqs. ( |2.8| , p.9| , p.l3| ), which lead to G N ^ 
find the renormalization group (RG) equations 



and the 



Z (N+N)/2 wethen 



N + N 
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Dp stands for the renormalization group differential operator 

T>» = fidp + (3 u d u + (t(k t - 2) + v 2 K VT ^d T + v(k v - l)d v 
Here we have introduced the Gell-Mann-Low functions 

A, du 



TK T + V K t 



din fi 

Or 



-e + 37 - 7„ )u 



<91n/i 

dv 



T 7-7, 



VK V 



u 7 - 7„ 



(3.2) 
(3.3) 

(3.4) 
(3.5) 
(3.6) 



and the Wilson- functions 7... = din Z.../dln fi\o. 

The RG equations can be solved in terms of a single flow parameter I using the charac- 
teristics 



/^(0=rJ(0(« o (fi(0)-i) 



it 



■u(l) = u 



Z^f (0 = (f Mu(0) - 2) + ^ T («(0) 
With help of these flow equations we recast Eq. ( |3.2D as 



f(l)=T 



G JV7 .({x},r(0^(/),w(/)},//i) = 



(3.7) 
(3.8) 
(3.9) 

(3.10) 



Equations (|3.7| - |3.9|) describe how the parameters transform if we change the momentum 
scale /t according to \i — > /i(/) = //i. Being interested in the infrared (IR) behavior of the 
theory, we must study the limit I — > 0. According to Eq. (|3.7| ) we expect that in this IR 
limit the coupling constant u(l) flows to a stable fixed point with f3 w = 0. At the fixed 
point it is legitimate to diagonalize the part of the RG differential operator containing the 
relevant parameters r and v. Introducing a new parameter a instead of r, we find 



VK v *d v 



(3.11) 



where 



The Ki* stand for the Gell-Mann-Low functions taken at the fixed point it*. 

Using dimensional analysis in conjunction with the flow equations, we readily find the 
asymptotic behavior of the connected Green's functions for I — > 0. Neglecting nonuniversal 
scale factors we get 

= ^lf N+ ^ d ^in( N ^l 2 G , (3.13) 

where the Fisher exponent r\ is defined by 77 = 7(11*). We define the remaining exponents by 

"=TZ?-. P = - d ~^ ±1 . (3-14) 



1 n d-2 + T] 

Vip = , Pip = vip - , (3.15) 

2 — K 7 



and the crossover exponent <fi by 

2 - 



1 



vjv 1P = p/frp . (3.16) 



The appropriate choice of the flow parameter I in the case |cr|<C|f|isi~|f| v . In this case 
the Green's functions scale as 

G^({x}) = \v\ iN+ ^F N ^({x/£},a/ \v\+) (3.17) 

with a correlation length £ ~ n~ l \v \ v . Otherwise, in the case |i>| <C lo") 1 ^, we choose 
/ ~ \a\ v and arrive at isotropic percolation scaling 

G^({x}) = Wl^^F^/M^/l^*) , (3.18) 
with a correlation length £/p ~ |cr|" /p . 



IV. RESULTS FOR THE CROSSOVER EXPONENT FROM THE e-EXPANSION 

In this section we derive the scaling indices. Because the e-expansions of the usual 
percolation exponents are well known, we concentrate on the new crossover exponent (p. For 
this purpose we need the Gell-Mann-Low functions (3 U , 7, and Ki explicitly. From Eq. ( p.4|) 
we know that f3 u begins with the zero-loop term — eu, the higher order terms are determined 
by the Wilson functions. These functions, the logarithmic derivatives of the Z-factors, are 
given by 7... = /id M | lnZ = (3 u d u \aZ r In minimal renormalization the Z-factors have 

a pure Laurent expansion with respect to e: Z = 1 + Y^/e + Y^/e 2 H . It follows 

then recursively in the loop expansion that the Wilson functions have also a pure Laurent 
expansion and, because they are finite for e — > 0, this expansion reduces to the constant 
term, i.e., all e-poles have to be compensated by the logarithmic derivative. Thus, we get 
the Wilson functions simply from the formula 7., = —ud u Y^ l \ Now it is easy to find these 
functions from the Z-factors Eqs. ( p.l0f]2~T2] . [2.2C| , |2.21| ). The results are 
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59 47 
lu = -Au + —u 2 + 0(u 3 ) , Tt — u + ——it 2 + 0(w 3 ) , (4.1) 
b 24 

7« = + + 0(« 3 ) , Jvr = ~ + + 0(U 3 ) , (4.2) 

3 21b 3 3b 

7 = -^ + |^ 2 + 0(« 3 ), (4.3) 

from which the Gell-Mann-Low functions follow as 

/5n=(-e + y-^ 2 + 0( M 3 )) M , . T = ^-^ 2 + 0( M 3 ), (4.4) 

= 1 ~ IT + > k vt = ^- ^u 2 + 0(u 3 ) . (4.5) 

b b 3 3b 

From (3 USf = the stable fixed point value it* = |e + J^£ 2 + 0(e 3 ) is readily obtained. We 
finally derive the following e-expansion of the crossover exponent: 

= 2-i + ^ + O(, 3 ). (4.6) 

Crudely evaluating this expansion for small spatial dimensions, i.e., for e = 3 or 4, leads 
inevitably to poor quantitative predictions. Therefore, we improve the e-expansion via a 
rational approximation, which takes into account that for d — 1 the correlation length 
exponents are, trivially, always equal to one. In addition, we make the hypothesis that the 
e-expansion can be extended up to e = 5. These considerations lead then to the interpolation 
formula 

* = 1 + V-V\ 1 + to + Ww)- (47) 

From Eq. ( |4.7|) we get values for the crossover and the order parameter exponent for d = 2 
to d — 6. These values are summarized in the following table: 



d = 


2 


3 


4 


5 


b 


= 


1.29 


1.52 


1.70 


1.8b 


2.00 


p = 


0.18 


0.b2 


1.09 


1.5b 


2.00 



Of course, the calculated values of <p depend slightly on the interpolation procedure, i.e., 
different rational approximands may be used to incorporate <p{d = 1) = 1. We learn from 
this numerical sorceries that the displayed numbers may have a failure of roughly ±0.05. 

Now we come back to the question which of the order parameter exponents are seen in 
simulations by Inui et al. 0. From the scaling properties of the Green's functions G N fr({x}) 
we have learned that this depends on the relative behavior of the relevant parameters \a\ 
and \v\^ as functions of the microscopic probabilities for the conducting elements. Writing 
the deviations from the multicritical line as 5p = p — p c (r), 5p± = p± — r, we have the 
expansions 

a = a 1 (5p + + 5pJ) + a 2 5p + a 3 (5p + + 5pJ) 2 + a 4 5p + 5p^ + a 5 5p 2 + ... , (4.8) 
v = (5 P+ - 5pJ)(b! + b 2 5p + b 3 {5 P+ + , (4.9) 
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where the coefficients ai, h depend on the microscopic model. In the simulations of the 
pure two-dimensional random diode system Inui et al. set Sp + = —8p_ = 5r and 5p = 0. 
It follows that er ~ Sr 2 and v ~ 5r, and consequently \a\ / \v\ ~ |<5r| 2 ~^ <C 1 because <p < 2. 
Thus the exponent (5 is found as the simulations clearly state. In their extended model Inui 
et al. set the variations of the probabilities to 5p + = 5p- = 5r and bp = 0. In this case we 
have a ~ 5r and v ~ 5r, and consequently \v\ / |o"| ~ fir] 1 - 1 '* < 1 because 1/0 < 1. 
Now the isotropic percolation exponent /3jp = 5/36 is measured. 



V. THE CROSSOVER FROM ISOTROPIC TO DIRECTED PERCOLATION 

A. Preliminaries 

The crossover theory of Frey, Tauber, and Schwabl 0, henceforth called FTS, starts 
with the quasistatic Hamiltonian 7i stated in Eq. ( [2.7|) . FTS use the modification of the 
longitudinal length scale by the parameter c which we have discussed in Sec. II. We have 
seen that c is a redundant variable near the isotropic fixed point and hence dropped it from 
the onset. FTS renormalize c also. They find that the anomalous dimension of c is zero at 
the isotropic point. Moreover, they calculate a nonvanishing anomalous dimension of c at 
the directed fixed point. The variable c (or better 1/c 2 , which couples to the composed field 
<y5V 2 <£>), being redundant at the isotropic fixed point, changes over to an irrelevant variable at 
the directed fixed point. Thus, in order to renormalize it, one has to include all the irrelevant 
operators of equal naive dimension and symmetry which mix under renormalization. This 
was overlooked in the work of FTS, thus the calculated anomalous dimension of 1/c 2 is 
meaningless. 

The technical problem of the crossover from IP to DP is, as FTS stated, that the two 
fixed points, which are to connect by the renormalization flow, have different upper critical 
dimensions, 6 and 5, respectively. Thus, one has from the outset the problem to renormalize 
the theory for general dimensions below 6 and cannot use the e-expansion, which by-passes 
the problem that the perturbation expansion is ill-defined if one uses critical massless prop- 
agators, leading to infrared (IR) divergencies below the upper critical dimension. However, 
the renormalization can be accomplished in a massive theory, which avoids the IR diver- 
gencies, by using normalization conditions for the vertex functions. Asymptotic scaling 
properties follow from the inhomogeneous Callan-Symanzik equation. 

Below the upper critical dimension, the theory is super-renormalizable. Only the vertex 
function Tn is UV divergent and can be renormalized by a mass shift r — ► r + 5t, where 
5t absorbs the UV divergencies. If the theory is regularized dimensionally, these UV di- 
vergencies manifest themselves as poles at e = Si = 2/1, I = 1,2,... p5|-{2~7]1, where, at the 
corresponding spatial dimensions di = 6 — si, only perturbational contributions to T^i of 
loop-order smaller then / are superficially UV divergent. These poles can be eliminated, 
as long as v is finite, by a shift 5t = g i ^ £ Ai(e), where Ai is meromorphic in e. After the 
mass shift, the theory is UV-finite. Then, however, in those dimensions corresponding to 
the before mentioned poles, non-analytic logarithmic behavior with respect to the coupling 
constant arises. 

An alternative dimensional regularization formalism was presented some time ago by 



Schloms and Dohm 241, henceforth called SD. This formalism circumvents renormalization 
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conditions and resorts to the more convenient minimal renormalization but without using 
the ^-expansion. The key feature of that formalism is to use the inverse correlation length 
itself as the mass-parameter [^7] and to define the minimal renormalization at a suitable 
renormalization point, which introduces the usual external mass scale /i. The consequences 
are relatively simple renormalization factors and a homogeneous RG equation. Hence, the 
SD formalism fosters higher order calculations and resummation procedures. We will use 
the SD formalism in the following. FTS use a massive renormalization scheme that is a 
mixture of renormalization conditions and a mass renormalization at an external mass scale 
fi. It is unclear wether this scheme can be used consistently at higher loop-order. 



B. Crossover Renormalization 



We revisit the quasistatic Hamiltonian TC, Eqn. (|2.7|) , and introduce, in the case v ^ 0, 
new variables by 

x ± = y j X|| = n ■ x = 2vXt , (p = |2i>| -1 / 2 s , tp — \2v\~ 1 ' 2 s , g — |2t>| 1//2 g . (5.1) 
Then 7i appears as 

1 „o . / _o\ Xg 



n 



dtd d - l y{ s 



s '-^ 2 + A ( T - v2 ) + f( s - s 



(5.2) 



In the limit v — > oo we arrive at the dynamic functional which describes directed percola- 
tion (DP) [[28|,[2lJ. We conclude that the so far unrenormalized Green's function have the 
asymptotic property 



{2v) ^N),2 G ({x}; T) V} g) _ Q ({x±; Mh ^ - g) 



(5.3) 



where the G N ^ are the DP Green's functions. The minimal renormalizations of the DP 
theory in dimensional regularization are 

s — > s = ZJ s , s — > s = ZJ 



r s , A — ► A — Z v ZX 



, . g 2 -> <? 2 = G E 1 Z V 1 Z 2 Z u u^ £ 



T — > f = Z Z T Tfi + T c 

where e = e - 1 = 5 - d = 4 - d, G s = T(l + e/2)/(47r) d/2 , and known 

■ 2 '13 31 35, 4* 



as 



Z = 1 



it 



tt u z /7 9, 4\ ^ ., 

1 + — + — (--3 + -In- )+0(u 3 ) 



s 128e V e 
\: ' 32e\e " ' 2 ""' 3 

{--- 

2e 2e \e 16 



+ — ln-l +0(u 3 ) 
4 2 3/ 



(5.4) 
(5.5) 



(5.6) 
(5.7) 
(5.8) 
(5.9) 



Note that the DP-Hamiltonian Ti, Eq. (|5.2|) , is non-renormalizable for d > 5. Thus, to study 
this limit we are restricted to spatial dimensions d < 5. 



- , u u 2 (\ 
Z T = 1 + — + 



0(u 



, 2m 7u 2 /! x A ^/ 3n 
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We expect that the theory can be rendered finite for all values of v, including infinity, 



by interpolating the renormalizations of Ji. Then it follows from the definitions (5.1) that 
the following asymptotic properties hold for v — > oo: 

Z...{u,v,e) -»• Zju,e) , (5.10) 

with u = u/2v. Of course, as long as v is finite the soft renormalizations of 7i can be used, but 
for v — > oo, logarithmic divergencies ~ \nv arise for e = e — 1 = 0. These must be included 
in the renormalization factors to make the theory finite for all v , especially in the DP limit. 
Note that in this limit v plays the role of a cutoff. In the following we demonstrate that 
the renormalization procedure suggested by SD, i.e., an extended minimal renormalization 
scheme without using the e-expansion, is feasible and that it allows to extract the crossover 
behavior of the Green's functions. 

We consider the dimensional regularized bare vertex functions N ({q};r,v,g;d) as 
functions of the bare parameters f,v,g, and the momenta {q} at spatial dimension d = 6—e, 
e > 0. The critical point f c is determined by 

f 1)1 ({0};f c ,v^;d) = 0, (5.11) 

which provides an implicit definition of 

r c = r c (v, g; d) = g"<S(e, v/g 2 ^) . (5.12) 



The generalized Symanzik function p5| has pole singularities at s — 2/1, I — 1, 2, as long 
as x = v/g 2 ^ £ is finite. An expansion yields S(s,x) = S (s) + Y(e)x 2 + 0{x A ) with a finite 
Y(e) as long as e > 0. Taking the limit behavior of the vertex functions for x ^> 1 into 
account, we find for e > 1, the asymptotic behavior S(e,x) — > x 2 ~ £ P(x~ £ ; e) + x 2 ^ 1 ^^ S(e), 
where P denotes an analytic function of x~ e = g 2 /v e . For e > 1 , the pole-singularities of 
P(y;e) and S(e) are found at e — 1 = e = 2/1. These pole-singularities combine to yield 
logarithmic divergencies in x instead the poles in e. 

Since the theory is super- renormalizable for e > 0, we know that all the functions 

f 5v )Ar ({q}; ^ v, g; d) = f ^ j7V ({q}; f' + f c (v, g; d), v, g; d) (5.13) 

are finite below six dimensions. However, because of the non-analyticity of f c with respect to 
g, they are no more expandable in the coupling constant g. This can be fixed by introducing 
the inverse transversal correlation length as the new mass. It is defined according to 

c-2 2 ^'iA{0};T',v,g;d) 

= m = 3 . (5-14) 

d q lT' hl ({q};T',v,g;d) 

The function m(f ', v, g; d) has poles only at d — 6. It can be inverted to define f' as a 
function of m, 

t = r(m,v,g;d) , (5.15) 
with r(0, v,g; d) = 0. The function r can be substituted into r'~ yielding the bare vertex 



functions Y'~ in terms of the mass m: 

N.N 



12 



r'ft, N ({Qhr(rn,v,g;d),v,g;d) = f '~ N ({q};m,v,g;d) . (5.16) 

In the following we abbreviate f ~ ^ by for notational simplicity. The T~ N are now 
free of dimensional singularities below d — 6, and have an expansion in integer powers of g. 
In the limit v —>■ oo, for d < 5, we deduce from Eq. (|5.3|) the asymptotic property 



t^ N ({q},m,v,g;d) - ( 2 £)(^- 2 )/ 2 f ^({q±, A" 1 ^}, m, |; d) , (5.17) 

with o> = 2\vq\\ and the unrenormalized DP vertex functions N . 

Henceforth we use for simplicity the notation q symbolically for all the momenta {q}. It 
is convenient to write 

rft jN (q,m,v,g;d) = rn»- N F^^ N (q/m,v/m, g/m e/2 ;d) (5.18) 

with dimensionless functions F% N and 5~ — d — (d — 2)(N + N)/2. Expanding F^ N in a 
power series gives 

oo . 

F~ jN (x,y,z;d) = z° £ f^M(x) n (^/e) , (5.19) 

l,n=0 

where a = if (N + iV)/2 is an integer and a = 1 if not. The exponent I denotes the loop 



f (i,n), 
7V,Af 1 

in y and can be expanded to 



order. The functions f^fZiV'i d) are finite for d < 6 and finite y. Moreover, they are analytic 



oc 



k=0 



The contributions of order z l to the scaled vertex functions F^ N ( |5.19 ) exhibit infrared 



(IR) divergencies for m — > 0. SD showed at the instance of the $ 4 -theory, that this problem 
can be treated by minimal renormalizations 

% = Z 1/2 (p , = Z 1/2 <p , (5.21) 
g 2 = A(d)G^Z~ 3 Z u ufj, £ , v -> v = Z~ x Z v v^i , (5.22) 

where the soft renormalization factors Z__ (.u; d) absorb just the poles at d = 6, followed by 
the application of the renormalization group. The UV e-poles and the IR divergencies are 
then summed up yielding the correct critical scaling behavior reminiscent of the intimate 
relation between both. Since the poles of the coefficients /-'^(Q — e)/e l do not depend on 
any of the parameters, they are identical to the usual soft minimal renormalizations and can 
be calculated by setting m = fi or in the massless theory with e — >• 0. Following SD, we have 
introduced in Eq. ( |5.22| ) a further amplitude A(d) with A(0) = 1. This dimension dependent 
function can be conveniently defined and is extremely useful for the practical calculation of 
scaling functions. For simplicity we set A(d) = 1 in the remaining part of this article. Note 
that from its definition (|5.14| ), the mass m does not need a multiplicative renormalization 
and that m = rh. 
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Now we consider the DP vertex functions 



N.N 



(g ± ,A 1 u,m,g\d) =m n ' n Ft, (q ± /m, us/km 2 , g/nrf/ 2 \ d) 



with <5„ 

d < 5 



(d + 1) — (d — l)(iV + N) /2 . The power series expansion of the F 



N.N 



N.N 



(x±,x\\,z;d) = z u f 



Z,nj_,ni|=0 



* {l ' n ^\d){x L )^{x^(z 2 /e) 1 



N.N 



(5.23) 
reads for 

(5.24) 



with e = 5 — d = 4 — d. By the same argumentation as above, one can show that the 
e-poles are cancelled and the IR divergencies are summed up by the renormalization group 
equation based on the minimal renormalizations Z, Z v , and Z u , Eqs. (5TR5T7,ST^), where 
the renormalization factors Z...{u, d) now absorb just the e -poles at d = 5. As above, they 
are identical to the usual soft minimal renormalizations. They can be calculated at the 
renormalization point m = \x , {q = 0} or in the massless theory with e —>■ 0, because the 

functions f~ n ^' n ^ (J) do not depend on any of the parameters. 

Combining the asymptotic properties of the vertex functions for y oo, Eq. ( |5.17| ), with 
the various expansions, we find for d < 5 



f { ~' n) (vd) 

N N 



(N+N-2-a)/2+nn-l 



N,N 



(d). 



(5.25) 



In case of taking the limit d — *■ 5 before y — > oo, the £T _i -poles are replaced by a poly- 
nomial p( l ' n± > n \0(lny) of order I. Therefore, it should be possible to find further finite 
renormalizations C...( u > v \ d) which interpolate between the two minimal renormalizations at 
the renormalization point m = /i, 



Z...(u; d) — > v] d)Z„_(u; d) = Z___(u, v; d) 



(5.26) 



where the Z...(u, v; d) tend to the corresponding Z___{u; d) with u — u/2v in the limit v — > oo. 

Following SD, the remaining problem is to determine the renormalization of the relation 
between the mass m and the "temperature" f = f(m,v, g; d) = m 2 T(v/m, g/m e ; d). The 
dimensionless function T has an expansion 



f{y,z;d) = l + Y,ti{y;d)(z 2 /e) 
i=i 



(5.27) 



but the ti(y;d) show UV singularities below d = 6. Otherwise, the use of the variable 
f = tt c = r(m,v, g;d) = m 2 T(v/m, g/m e ;d) — g^ £ S(e,v/g 2 ^ e ) instead of f eliminates all 
UV singularities below the upper critical dimension. Then, however, the function r is not 
more expandable in g. Following SD, we shall therefore consider the derivative dr/dm 2 \ i ^, 
which is not only expandable but also free of singularities. We define the function 



P(v/m, g/m £ ; d) 
which has the expansion 



Or 



dm 2 



d{m 2 T) 



v,9 



dm 2 



Or 



",9 



dm 2 



(5.28) 



v.y 
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P{y,z;d) = l + Y / Pi(y;d)(z 2 /e) 
1=1 



(5.29) 



The functions pi(y; d) are finite for d < 6, and P is minimally renormalized by 

P = Z~ X Z T P . (5.30) 

The same argumentation as above shows that it is possible to find a Z T (u,v;d), which 
interpolates between the respective minimal renormalizations of isotropic and directed per- 
colation. Note that <9/<9m 2 |^ = d/dm 2 \ vg . Thus, the derivative with respect to the mass 
m commutes with the multiplication with Z-factors. 



C. One-Loop Crossover Calculation 

In this subsection we derive a minimal crossover renormalization, the corresponding 
renormalization group equation, and the flow equations, which show the crossover from 
isotropic to directed percolation. From the self-energy to order g 2 , we get the unrenormalized 
(for notational simplicity we drop the overcirc) vertex function I^i expanded to second order 
in the momenta: 

Ti,i = r 1 7- r — K^v/y/T) +2iv-q 1 K^ 2 (v 



(2-e)e 1 v V "V 2 £ 

Vfl - ^^(^) - (v ■ q) 2 ^^»(^) • (5.31) 



The inverse transversal correlation length is given by 

i _ a^/fW (t y^ 



>» 2 = >. G J;-iw»;:v^ ■ (5.32) 



From Eq. ( |5.32|) we find easily the unrenormalized perturbational temperature as 

r = m 2 (l + ^E^Q* - K^v/m))) . (5.33) 



Taking the derivative with respect to m 2 while holding g and t> constant and renormal- 
ization by multiplication with the factor ZZ' 1 yields via the renormalization scheme the 
function P as 

P = ZZ; 1 + u(fi/m) e (^(iK^} 2 (vfi/m) - K^vfi/m)) - ^K^(vfi/m)^j . (5.34) 



With the formulas 



p 2 K£\p) = (l+p 2 )K a (p) - K a „ 2 (p) , (5.35) 
I + a)K a (p) = (2 + a)K a+2 (p) - (1 + p 2 )- 1 ^ 2 , (5.36) 
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where we have defined K a := Kf \ we get 



P = zz; 1 + ^(^(i + v r £/2 



8\3e 



+ 



£- 1 



(3 - v-z) 



e-3 



-2 



K £ (v)-(l + v 2 y £/2 



-^(i + ,V £/2 



K e (v)-(l + v 2 y £/2 
+ KJv)) 



(5.37) 



at the renormalization point m = \i. Note that the poles for e — 1 and e = 3 are only 
fictitious because K 1 (p) = (l+p 2 )^ 1 ^ 2 and K 3 (p) = (l+p 2 ) _3 / 2 (l + 2p 2 /3). Useful properties 
of the function K £ {p) are 



K £ (p) = l- 2 -^ P 2 + 0(p*) 

_ y/HX(l + g )/2)^ 1 1 ^ 2 _ £ 



2r(l+e/2) 



1 + e 



(5.38) 
(5.39) 



from which we find that the fictitious pole at e — 1 in Eq. (|5.37|) is replaced by a divergence 
~ lnt> , and that w = uK £ (v) — > u/2 for v — > oo. Thus, a minimal crossover renormalization 
is given by 



ZZ- 1 = 1 - — (1 + t; 2 )- £/2 - -4^ 
6e v ' 4(e-l) 

Defining the crossover function 



K £ (t;)-(l+t; 2 ) 



2\-e/2 



0(w 2 



C e (v) = !-((! + v 2 ) e/2 /f e ( V ))" 



we render the renormalized temperature derivative 



w 



P = l + - I 



1-CJv) . 5-e 



3-e 



+ 



1-e 



(5.40) 



(5.41) 



(5.42) 



finite for all d < 5 and all v. 

Now we consider the vertex function T^i, Eq. ( |5.31|) , for q 7^ 0. Inserting r, Eq. ( |5.33|) , 
we get after renormalization 



r lil = (m 2 + g 2 )^ 

+2i/iv • q^Z„ - 



4e 



(v . q) . »(>-w\ ? ) Wro) (543) 



The last term is finite for all d < 5 and all t> , even in the limit u — > oo, if one holds 
2v-q=c<j/A and u> = uK £ (v) constant. The reduction of K^} 2 { v ) with Eqs. (|5.35| , |5.36D 
leads to 



2 + v 



-2 r 



# e (u)-(l+<) 



2\-e/2 



i^(i„-[^ w -(i+ B r e/2 



3<l+« 



2N-e/2 



(5.44) 
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We see that the divergencies which have to be absorbed by renormalization arise from the 
first two terms of the right hand side of Eq. ( |5.44| ). Thus, we define 



u 



u 



z „. 1 + Ji (1 + „ r ./ 2 + _JL_ 



K £ (v)-(l + v 2 y £/2 
K e (v)-(l + v 2 y £/2 



+ 0(u 2 ) , 
+ 0(u 2 ) . 



(5.45) 
(5.46) 



Now we turn to the calculation of the one- loop contribution to the vertex function r 12 . 
To renormalize T 12 we need it for zero external momenta. We find the contribution 



W = -2g* [ G(k) 2 G(-k) = 2^£ (1) (0) = ^r^ 2 K £ . 2 (v/ 77) 
A or e 



which leads after renormalization to 



(T 1<2 ) 2 = G; 1 u f i £ [Z u 



Au(n/m) 



•K e - 2 {vn/m] 



Using once more the reduction formula ( |5.36| ), we finally find 



An, 

Z u = l + ^(l+v 



2\s/2 



+ 



Au 



K £ (v)-(l + v 2 ) 



2\s/2 



0{u 2 



Now, all the Z-factors ( |5.4Q]5.45] . p.46| , |5~4*9"| ) are of the form 



z l = i +u ^(i+v 2 r /2 + 



[e-l) 



2\s/2 



where the coefficients are given by 

a = 1/6 , 6 = 1/4, 
a v = 1/3 , b v = 1/2 



K £ (v)-(l+v 2 ) 



a T — 1 , b T — 1 , 
a,. = A , b u = A . 



+ 0(u 2 ) 



With help of the derivative of the i^-function, 

(1 + v 2 )- l - £/2 - K £ (v) = -K £ (v)C £ (v) 



dKJv) 



(5.47) 



(5.48) 



(5.49) 



dv 

and the Gell-Mann-Low functions 

(3 U = du/d\nfi\ Q = (—e + 37 - 7„)u = —eu + 0(u 2 ) , 
(3 V = du/d\an\ = (-1 + 7 - Jv)v = (-1 + 0{u))v 

we obtain the Wilson-functions 7i = <91nZj/<91n/i| = {(3 u d u + (3 v d v ) lnZj| easily as 



7i = -u( Oi(l + v 2 ) 1 6/2 + bi 



K £ (v)-(l+v 2 ) 



2\-l-e/2 



+ 0{u 2 



-wU(l-C e {v)) + biC e (v)\+0(u 2 ) . 



(5.50) 



(5.51) 
(5.52) 



(5.53) 



(5.54) 
(5.55) 



(5.56) 
(5.57) 
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Here we have used the new coupling constant w = uK £ {v) and the definition Q5.41D of the 
crossover function C £ (v) . We mention that w = u for v = and w = u/2 for v = oo, as well 
as C e (0) = and C £ (oo) = 1. C £ {y) crosses over monotonically between this two values. 
Especially in three dimensions it is a simple rational function of v 2 : 



C 3 (v) 



v 2 (5 + 2v 2 ) 
;i + v 2 ){3 + 2v 2 ) ' 



(5.58) 



The renormalization group equation 



N + N 



-1 



G 



(5.59) 



follows as usually from the independence of the unrenormalized Green's functions from the 
external mass scale \x. Here we have to note that the mass m is by definition ( 5.14 ) a function 
of the bare parameters only. Therefore, its Gell-Mann-Low function vanishes. Thus, we have 
the renormalization group differential operator 



V. = fidf, + (3 w d w + (3 v d v 



(5.60) 



To one-loop order we obtain 



-e-C e {v) + 



42 - 5C £ (v) - C e {vf 



12 



w \w 



, 2 + C £ (v) 

-1 H -^w \v 

12 



and 



7 



2 + C £ (v) 



12 



-w . 



The flow equations 



lj l w(l) = (3 w (w(l),v(l)) 



w(l) = w 



lj l v(l) = /3 v (w(l),v(l)) , v(l) = v 



(5.61) 
(5.62) 

(5.63) 

(5.64) 
(5.65) 



and their solutions show of course the (unstable) isotropic and the (stable) directed percola- 
tion fixed points for v * = and v* = oo, respectively, and the continuous crossover between 
both. 

Upon defining the flowing amplitude function X(l) by 



i|lnX(i)= 7 (tS(0,«(0) 



X(l) 



(5.66) 



the scaling form of the Green's functions follows from the renormalization group equa- 
tion ( |5.59|) as 
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K (N+N)/2 



" ( md ~ 2 ^(m/^ ' (5 ' 67) 



where the last form holds asymptotically for v(m/fi) ^> 1, and A/A(/) = lv(l)/v. 

We conclude this section with the determination of the parameter m as a function of the 
temperature r. For the renormalized temperature derivative we have the renormalization 
group equation 

T? M - k P = (5.68) 

with k = 7 — 7 r . Its solution is given by 

P{m/ fi;w,v) = Y{m/ fi)P(l;w(m/ fi),v(m/ /i)) , (5.69) 
where the amplitude function Y is determined by the differential equation 

lf^nY(l) = -k(w(1),v(1)) , Y(l) = 1 . (5.70) 



P comprises the main (exponential) contribution of the crossover-scaling of dr/dm 2 
The amplitude P(l;w,v) is given to one-loop order by Eq. ( p.42| ). The temperature r 
results then from the solution of the differential equation dr/dm 2 \ WVfl = P(raj fi;w,v). 
The integration of all the flow-equations has to be done numerically and leads, qualitatively, 
back to the numerical results of FTS 0. 



VI. EPILOGUE 

Using field theoretic methods, we have analyzed the connectivity behavior of random 
resistor-diode networks near the percolation critical point. We found that the introduction 
of positive and negative diodes oriented to a privileged direction in space with unequal 
probabilities leads to a crossover to the directed percolation problem, whereas a distribution 
of diodes with equal probabilities results only in elongated isotropic percolation clusters. In 
the latter case, a simple rescaling of the privileged direction maps the problem to isotropic 
percolation. A slightly different distribution of the diodes introduces a further relevant 
variable with a new scaling dimension. We have calculated this scaling dimension to 0(e 2 ) 
in an e-expansion around six dimensions. An interpolation resulting from this e-expansion 
and an exact value at one dimension leads to a formula that compares very well with a recent 
simulational result in two dimensions. It would be very interesting to perform simulations 
also in higher dimensions in order to compare these with our result. 

In Sec. V we have reconsidered the theory by Frey, Tauber, and Schwabl for the crossover 
from isotropic to directed percolation. Some shortcomings are corrected and it is demon- 
strated, how one can perform crossover calculations consistently by using a type of extended 
minimal renormalization. 
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APPENDIX: TWO-LOOP CALCULATION 



In this appendix we present briefly the main part of the two-loop calculation of the new 
renormalization constants Z v and Y VT . The dimensionaly regularized parameter integral 



I{a,b, c) 



P,q (a + p 2 )(6 + q 2 )(c+(p + q)2) 



25 



Qe 12 
^3 21 



_3— e i l3— e 



+ b 6 ' 6 + c 



3-e\ 



— 3abc 



~(J + t) (fl2 ~ £(6 + c) + 62 ~ £(a + c) + c2_£(6 + c)) ) ' 

introduced in Ref. [50f, plays a fundamental role in the calculation. Its derivatives 

(_l)^+™+n-3 d l+m+n ~ 3 I{a, b, c) 



L lmn 



(I - l)!(m - l)!(n - 1)! da^db^dc^ 1 



a=b=c=l 



are extensively used in the following. Particularly simple are the integrals 



1 



d n-l a 2-e/2 



P (l + p 2 ) n (n-l)!(4-e)(2-e)e da n ~ l 



(Al) 



(A2) 



(A3) 



a=0 



We start with the self-energy diagram displayed in Fig. 2(a). Its value is given by 



E (2a) (q) = / G(q - k)G(k) 2 G( P )G(k - p) 



1 



2 7k, P [f + (q-k) 2 ][f + k 2 ] 2 [f + ( p+ w) 2 ][r+(k-p 



,21 ' 



(A4) 



where we have shifted k — > k— zv, and defined r = r + v 2 and q = q + 2iv. We expand the 
last integral in q and v to second order and find 



S (2a) ( 



-d-4 



Af 



t( /ii3 + (V /i23 + t /i33 ) v2 



i\ ■ q 



d-A 



-h 



Af 
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d \ d 

In the same fashion we get for the diagram shown in Fig. 2(b) 

S (26) (q) = ~9 4 I G(q " k)G(q - p)G(k)G(p)G(k - p) 

Jk,D 



115 q 



(A5) 



/k, P 

-g~T \ 1 122 



i - - / 1 , 

-M9,9. + ( ; -<222 H T-'223 J v ~ 



d 



4f 



d-6 r 4f r 1 2 r_ 

~ /l23 + T /l24 + / 3 + d /l33jcr 



(A6) 
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Using Eqs. ( |A~2| , [A~3|) we obtain from Eqs. ( |A5|JA6|) the singular contributions to the two-loop 
vertex function after renormalization according to the scheme (^78, 2.9 , 2.13 ) as 



p(2— loop) 
1 1,1 



U 



/9 45e 



fl T 



5 25e A (/iv) z 
2 + ~2A 



11 



7e A q 2 , /23 13£x2^v-q 



72^ 6 



+ (l2 + 



12 144 



(A7) 



Using the renormalization constants to 0(u), Eqs. ( |2 . 1 0| - |2~T2"| , |2 . 2 0| , |2~2T| ) , we get from 



Eq. (|2.18|) , after the renormalization f x,i — > IYi = ZT^i, the one- loop self-energy to 0{ 



u 



as 



loop) 

1 1,1 



-e/2 



9 llexM 

-11 5e \ u \ q 2 
T ~ 12'£7~6 



/i 2 r 



1 



5e\ m 



12' £/ 

23 5e^ u\ 2ifiv • q 
6 



+ ( 1 + (T-n)7)l + ( 1 + (lT-^)7 

r (2- 
1 1,1 



12' e,/ 3 

By adding Y^ loop \ p( 2 ~'°°p) an( j ^ e renormalized zero-loop part, we get 



li 2 (Z T T + Y VT v 2 ) + 2iZ vt iv ■ q+Zq 2 + T^ loop) + T^ loop> + O(s , u 



-\(2— loop) 



-0 „,3\ 



(A8) 



(A9) 



We see that the nonprimitive divergencies ~ In r of loop ^ and loop ^ cancel (as a 
check of a correct calculation) and find finally the renormalization constants cited in 
Eqs. (|2~lQ| , pl] , ^20i^2lD . 
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FIGURES 



aIP+-RJ 




FIG. 1. Generic phase diagram for p + +p_ = const < 1— p. The nomenclature is the following, 
p c : critical point, nonP: non-percolating phase, DP: phase of directed percolation, IP: phase of 
(elongated) isotropic percolation. 
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